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Abstract 

This  paper  analyzes  an  immersed  finite  element  method  based  on  piewise  linear  polynomials 
introduced  recently  for  solving  interface  problems.  As  an  important  feature,  this  method  allows 
the  interface  to  pass  through  triangles  so  that  partitions  with  simple  structures,  such  as  the 
Cartesian  partition,  may  be  employed  to  solve  a  problem  with  an  arbitrary  interface.  The  main 
objectives  here  are  to  demonstrate  the  approximation  capability  of  the  immersed  finite  element 
space,  and  to  derive  error  estimates  for  the  related  finite  element  solutions.  Numerical  results 
are  also  provided  to  illustrate  some  properties  of  this  method. 

1  Introduction 

In  [22],  two  (nonconforming  and  conforming)  immersed  finite  element  methods  based  on  Carte¬ 
sian  grids  are  developed  for  boundary  value  problems  of  elliptic  differential  equations  with  dis¬ 
continuities  in  the  coefficients;  singularities  in  the  source  terms;  and  arbitrary  interfaces  in  the 
solution  domains.  These  problems  are  referred  as  interface  problems.  In  this  paper,  we  investi¬ 
gate  the  nonconforming  immersed  finite  element  method,  especially  the  related  error  estimates, 
in  a  more  general  framework. 

The  model  interface  problem  of  concern  consists  of  an  elliptic  equation  of  the  form 

(1.1)  -V  ■  (0Vu)  +  ku  =  /, 

where  3  has  jumps  across  an  interface  T  in  the  solution  domain  f l  £  R2,  together  with  the 
Dirichlet  boundary  condition.  In  addition,  we  allow  the  source  term  /  to  have  a  Dirac  delta 
function  singularity  on  T  of  the  form 

(1.2)  /(x)  =  /c(x)  +  <2(X)  <j(x  —  X)dX, 

where  <2(X)  is  the  source  strength;  <5  is  the  Dirac-delta  function  which  is  defined  in  the  sense 
of  distribution.  Such  a  source  function  is  an  important  feature  of  Peskin’s  immersed  boundary 
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method  (IBM)  [27,  28]  that  has  been  used  for  many  problems  in  mathematical  biology  and 
computational  fluid  mechanics  [3,  11,  12,  13,  29]. 

From  equation  (1.1)  and  (1.2),  it  is  easy  to  obtained  the  jump  conditions 

(1.3)  [«]  =  0,  continuity  condition, 

(1.4)  [/3it„]  =  Q(s),  net  flux  across  the  interface, 

where  the  jump  is  defined  as  the  difference  of  the  limiting  values  from  each  side  of  the  interface 
assuming  the  solution  is  piecewisely  smooth  on  each  side;  un  is  the  normal  derivative  of  the 
solution.  Many  applications  involve  solving  interface  interface  problems  similar  to  the  one 
defined  by  (1.1),  (1.3),  and  (1.4),  for  example,  the  projection  method  for  solving  the  Navier- 
Stokes  equation  [1,  2,  8],  two  phase  flow  [30],  Hele-Shaw  flow  [14,  15],  and  many  others. 

It  is  well  known  (see  [5,  7]  and  the  references  therein)  that  the  standard  Galerkin  method 
with  linear  finite  elements  can  be  used  to  solve  such  elliptic  interface  problems.  However,  to 
achieve  the  optimal  0(h2)  accuracy  in  the  numerical  solutions,  triangles  are  required  to  be 
aligned  with  the  interface,  i.e.,  the  interface  is  allowed  to  pass  a  triangle  only  through  its 
vertices.  This  restriction  will  obviously  prevent  the  Galerkin  method  with  linear  finite  elements 
to  work  efficiently  for  those  applications  in  which  an  interface  problem  similar  to  the  one  defined 
by  (1.1)-(1.4)  has  to  be  solved  repeatedly,  each  time  with  a  different  interface  T,  because  the 
partition  has  to  be  formed  over  and  over  again  due  to  the  variation  (either  the  shape  or  the 
position)  of  the  interfaces. 

In  the  conventional  Galerkin  finite  element  method,  the  linear  polynomial  basis  functions  in 
each  triangle  are  independent  of  the  interface  conditions,  but  the  partition  is  formed  according 
to  the  interface  so  that  a  better  numerical  result  can  be  achieved.  The  immersed  finite  element 
method  employs  an  opposite  approach:  the  partition  in  the  immersed  finite  element  method  can 
be  independent  of  the  interface,  but  the  basis  functions  are  constructed  according  to  the  jump 
conditions  across  interface.  The  immerse  finite  element  method  therefore  has  some  advantages 
in  solving  interface  problems.  First,  partitions  with  simple  structure,  for  example,  a  Cartesian 
partition,  can  be  used  even  if  the  interface  T  is  arbitrary.  The  only  constraint  on  a  partition  is 
that  it  must  be  fine  enough  to  satisfactorily  represent  the  interface.  Secondly,  it  is  possible  to  use 
one  or  just  a  few  partitions  to  solve  a  sequence  of  interface  problems  so  long  as  the  interfaces  do 
not  vary  too  much.  Also,  the  immersed  finite  element  method  may  be  used  together  with  other 
Cartesian  grid  methods  such  as  the  level  set  method  [26]  and  those  implemented  in  Clawpack 

[17]- 

Since  the  immersed  finite  element  method  has  no  restriction  on  the  partition,  triangles  in  the 
partition  can  be  separated  into  two  groups;  the  first  group  consists  of  interface  triangles,  and  the 
second  contains  noninterface  triangles.  A  triangle  is  an  interface  triangle  if  it  is  separated  into 
two  subsets  by  the  interface;  otherwise,  it  is  a  non-interface  triangle.  Obviously,  we  can  just  use 
the  usual  linear  basis  functions  to  construct  approximation  in  a  non-interface  triangle.  However, 
in  an  interface  triangle,  we  introduce  basis  functions  with  linear  polynomials  piecewisely  defined 
in  its  two  subsets  formed  by  the  interface  such  that  these  basis  functions  can  satisfy  the  jump 
conditions  (either  exactly  or  approximately)  on  the  interface  and  retain  specified  values  at  the 
vertices  of  this  interface  triangle.  The  idea  here  is  similar  to  that  used  for  the  Hsieh-Clough- 
Tocher  macro  C 1  element  [4]  where  each  basis  function  consists  of  three  cubic  polynomials  on 
the  sub-triangles  formed  by  connecting  the  vertices  and  the  center  of  gravity  so  that  the  required 
continuity  can  be  satisfied.  We  refer  the  readers  to  [10,  14,  18,  19,  20,  21,  23,  24,  25,  31]  for 
more  background  materials  about  immersed  interface  and  immersed  finite  element  methods,  as 
well  as  their  applications.  Our  main  effort  here  is  devoted  to  error  estimation  for  this  method. 

This  paper  is  organized  as  follows.  Section  2  contains  some  preliminary  materials  such  as  the 
specific  interface  problem  to  be  considered,  basic  notations,  an  immersed  finite  element  space 
for  the  interface  problem,  and  some  of  its  basic  properties.  Section  3  is  devoted  to  an  immersed 
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finite  element  method  for  the  interface  problem  and  the  related  optimal  error  estimates.  Section 
4  presents  a  modified  immersed  finite  element  method  whose  implementation  is  easier  from  the 
point  of  view  of  numerical  integration.  Section  5  contains  numerical  results  illustrating  features 
of  these  immersed  finite  element  methods. 


2  The  immersed  finite  element  (IFE)  space 

2.1  Preliminaries 

To  be  specific  but  without  loss  of  generality,  we  consider  the  immersed  finite  element  (IFE) 
method  for  the  following  boundary  value  problem: 

(2.5)  -V-(/3 Vu)=/,  (*,!/)  en, 

(2.6)  u\ an=g, 


together  with  the  jump  conditions  on  the  interface  T: 


(2.7)  Mir  =  0, 

(2.8)  [I3un]\r  =  0. 


Here,  see  the  sketch  in  Figure  1,  0  C  R2  is  a  convex  polygonal  domain,  the  interface  T  is  a 
smooth  curve  separating  fi  into  two  domains  fl_,  fl+  such  that  fl  =  fl_  U  Q+,  and  the  coefficient 
j3  (x,y)  is  a  piecewise  constant  function  defined  by 


(x,y)  e  n  , 
(x,y)  G  fi+. 


We  also  assume  that  the  boundition  function  g  is  zero,  and  the  most  of  the  results  obtained  in 
this  paper  can  be  extended  by  the  usual  way  to  the  cases  in  which  g  yf  0. 


Figure  1:  A  sketch  of  the  domain  for  the  interface  problem. 


We  let 

(2.9)  Hm(n)  =  {u  |  u|n-  £Hm{ns),  s  =  -,+}, 
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and  we  define  the  norm  of  to  be 


(2-10)  —  llUllm,£2-  +  IMIm,fi+  ■ 

Semi- norms  of  can  be  defined  accordingly  by 

(2-11)  Mm, £2  =  M  m,n-  +  Mm, £2+  ■ 

For  any  subset  T  of  fl,  we  let 

Ts=Tnffs,  «  =  -,+. 

For  any  function  f(x,y )  defined  in  T  C  fl,  we  can  restrict  it  to  Ts,s  =  — ,  +  to  obtain  two 
functions  as 

/*(*> y )  =  f(x, y),  if  (x, y)  e  Ts,  «  =  —,+. 

We  use  HF  to  denote  the  line  segment  between  two  points  D,E  £  fl.  For  any  curve  T,  we  use 
nr  to  denote  its  unit  normal  vector  pointing  to  a  particular  side  of  T.  Also,  for  any  measurable 
subset  A  of  fl,  we  use  |A|  to  denote  its  measure. 

2.2  Finite  element  functions  in  interface  triangles 

Let  Th  be  a  regular  partition  of  fl  with  a  step  size  h.  Without  loss  of  generality,  we  assume  that 
the  triangles  in  this  partition  have  the  following  features: 

{Hi):  If  T  meets  one  edge  of  a  triangle  at  more  than  two  points,  then  this  edge  is  part  of  T. 

{H2):  If  T  meets  a  triangle  at  two  points,  then  these  two  points  must  be  on  different  edges  of 
this  triangle. 

For  a  typical  triangle  T  £  Th,  we  use  A  —  {xi,y2),B  =  {. x2,y2),C  =  (£3,2/3)  to  denote  its 
vertices,  and  use  D  =  ( xn,yD )  and  E  =  {xe,  Ve)  to  denote  its  interface  points  if  T  is  an 
interface  triangle,  see  the  sketch  in  Figure  2. 


Figure  2:  A  typical  interface  triangle  A  ABC.  The  curve  between  D  and  E  is  part  of  the  interface 

r. 
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We  now  consider  the  finite  element  functions  in  an  interface  triangle.  We  follow  an  idea 
similar  to  that  for  the  Hsieh-Clough-Tocher  macro  Cl  element  [4]  in  which  piecewise  polyno¬ 
mials  are  used  in  a  triangle  to  maintain  certain  desirable  features.  For  our  interface  problem, 
we  obviously  would  like  the  finite  element  functions  to  satisfy  the  jump  conditions  across  the 
interface.  Since  the  interface  T  separates  an  interface  triangle  T  into  two  subsets  T~  and  T+,  we 
naturally  can  try  to  form  a  finite  element  function  by  two  first  degree  polynomials  defined  in  T~ 
and  T+,  respectively.  Note  that  each  polynomial  of  degree  one  has  three  freedoms  (coefficients). 
The  values  of  the  finite  element  function  at  the  vertices  of  T  provides  three  restrictions.  The 
normal  derivative  jump  condition  provides  another.  Then  we  can  have  two  more  restrictions  by 
requiring  the  continuity  of  the  finite  element  function  at  interface  points  D  and  E.  Intuitively, 
these  six  conditions  can  yield  the  desired  piecewise  linear  polynomial  in  an  interface  triangle. 
This  leads  us  to  consider  functions  defined  as  follows: 


(2.12) 


4>{x,  y)  =  { 


4>  {x,y)  =  axx  +  biy  +  cx, 

<t>+  (, x ,  y)  =  a2x  +  b2y  +  c2, 

4>-{d)  =  4>+{d),  cp~(E)  =  4>+(E), 

rvr-w)-nM=o, 


(x,y)  e  T  , 

(*,!/),€  T+, 


where  is  the  unit  vector  perpendicular  to  the  line  DE. 

Lemma  2.1  The  function  <j>(x,  y)  defined  by  (2.12)  in  an  interface  triangle  T  is  uniquely  decided 
by  its  values  at  the  three  vertices  ofT. 

Proof.  We  need  to  prove  this  only  in  the  reference  triangle  T  whose  vertices  are  Ai  =  (0, 0)T,  A2  = 
(1,0)T  and  A3  =  (0, 1)T.  Without  loosing  generality,  we  assume  that  <j>{x,y)  takes  a  value  di  at 
the  vertex  Aj  for  i  =  1,2,3,  respectively,  and  the  interface  T  meet  this  triangle  at  E  =  (a,  0)T 
and  D  —  (0,  b)T  such  that  A\  €  T~  with  0  <  a  <  1, 0  <  6  <  1,  and  ab  <  1.  The  values  specified 
at  the  vertices  force  <f>(x,  y)  to  have  the  following  form: 


,(  a  _  f  {x,y)  =a1x  +  b1y  +  d1,  {x,y)  ef  , 

\  4>+(x,y)  =  d3{l  -  x)  +  d2(l  -  y)  +c3(l  -  x  -  y),  {. x,y),£T+ , 

with  ai, bi  and  C3  to  be  determined.  The  continuity  condition  at  points  E  and  D  leads  to 


(2.13) 


J  a,ia  +  di  = 

(  bib  +  d\  = 


—  ^3(1  —  a)  +  d2  +  03(1  —  a), 
ds  +  d2(l  —  b)  +  03(1  —  b). 


The  last  equation  in  (2.12)  requires 


(2.14) 


^(ttirii  +  b\n2 )  —  —  [( d's  +  C3)n±  +  (d2  +  03)712] , 


where  n  =  (ni,n2)T  the  unit  vector  perpendicular  to  the  line  DE  .  Therefore,  the  coefficients 
ai,  bi  and  C3  of  f>  satisfies  the  following  linear  system: 

e?3  (1  —  a)  +  d2  —  di 
ds  +  d2(l  —  b)  —  d\ 

-d^rii  -  d2n2 


r 

/?+' 


(  a 

0 

a  —  1  ^ 

1  (  Ql 

0 

b 

6-1 

bl 

\  Rni 

Rn2 

711+712  ) 

'{cs 

with 

(2.15) 

R  = 
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A  simple  calculation  shows  that  the  determinant  of  the  matrix  in  this  system  is 
ab  jni  +  n2  +  R  [(1  —  a)  ^  +  (1  —  b)  ^ j  j  , 

which  can  not  be  zero  since  n  =  (ni,ri2)T  is  not  a  zero  vector  and  rii .  and  ri2  do  not  have 
different  signs.  Thus,  the  above  linear  system  must  have  a  unique  solution,  and  the  function 
c p(x ,  y)  here  is  uniquely  determined  by  its  values  at  the  vertices  of  an  interface  triangle  T. 


We  now  let  Sh(T)  be  the  linear  space  of  all  the  functions  defined  by  (2.12),  and  call  it  the 
immersed  finite  element  space  in  an  interface  triangle  T.  The  lemma  above  shows  that  the 
dimension  of  Sh{T)  is  three.  Furthermore,  we  notice: 

•  The  proof  of  Lemma  2.1  itself  provides  a  way  to  construct  the  nodal  basis  functions  in  an 
interface  triangle. 

•  From  the  proof  we  can  see  that  cj)~(x,y)  =  cj)+(x,y)  when  R  =  1,  i.e.,  when  the  coefficient 
does  not  have  the  jump,  the  functions  of  S)j(T)  become  the  usual  linear  polynomials.  In 
this  case,  S/,(T )  reduces  to  the  standard  linear  finite  element  space. 

•  When  rnT  is  a  straight  line,  the  function  cf>(x , y)  defined  by  (2.12)  is  continuous  in  T  and 
therefore  is  in  U1(T). 

2.3  Interpolation  errors 

We  now  discuss  the  approximation  capability  of  the  space  S(T)  when  T  is  an  interface  triangle, 
and  we  will  follow  the  usual  framework  for  the  theory  of  affine  families.  Consider  a  set  J(T)  of 
functions  such  that  every  u  £  J(T)  satisfies: 

(  u\Ti  =«'  e  H2(Tl ),  i  =  -,  +, 

(2.16)  u~(D)  =u+(D),  u~{E)  =u+{E), 

l  JrnT^~^u~  ~  P+^u+) '  nrds  =  0. 

The  set  J(T)  is  a  linear  space  in  the  usual  sense,  and  for  each  non-negative  integer  m  <  2,  we 
endow  this  space  with  a  norm  ||-||m  T  defined  in  the  same  way  as  (2.10).  Semi-norms  of  J(T) 
are  defined  similarly. 

Lemma  2.2  For  an  interface  triangle  T,  the  space  Sh(T)  is  a  subspace  of  J(T). 

Proof.  For  any  cj)  £  Sh{T),  it  is  obvious  that  (ps  £  H2(TS),  «  =  —,+.  Also,  because  ^  is  a 

piecewise  linear  polynomial  satisfying  (2.12),  Green’s  formula  leads  to 

[  (/T V<T  -  /?+V^»+)  •  nr ds  =  -  /  -  /?+V^+)  •  n^ds  =  0. 

irnT  Jde 

Thus  (j>  £  J(T)  and  Sh(T)  is  a  subspace  of  J(T). 


Remark:  As  a  consequence  of  this  lemma,  any  immersed  finite  element  function  (f>  £  Sh{T)  can 
satisfy  the  derivative  jump  condition  exactly  in  a  weak  sense  as  follows: 

[  {/3~V(f>~  -  /3+V<f>+)  ■  nrds  =  0. 

JrnT 

Now,  for  any  u  £  J(T)  we  let 


(2.17) 


12,  T 


|u|2  T  +  |u(A)|  +  \u(B)\  +  |u(C)| . 
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Lemma  2.3  |||  •  |||2,t  is  a  norm  in  the  space  J(T),  and  this  norm  is  equivalent  to  ||-||2  T. 

Proof.  We  mainly  need  to  show  the  equivalence.  By  the  embedding  theorem  of  Sobolev  spaces, 
it  is  easy  to  see  that  for  any  u  £  J(T),  there  exists  a  constant  C±  >  0  such  that 

(2.18)  |||u|||2,t  £  C±  1 1 1 1 ^  j’  ■ 

Then  we  would  like  to  show  that  there  exists  another  constant  Ci  >  0  such  that 


(2.19)  IMI2  T  —  ^2||M||2,T- 

Assume  that  this  is  not  true,  then  there  exists  a  sequence  {u&}  c  J(T)  such  that 

(2.20)  IMI2,t  =  i,  IIKIIl2,T<p 

and  we  can  further  assume  that  this  sequence  converges  in  1 1  ■  1 1  x  j,  norm  of  J(T).  Since 

||u;  —  Vm  ||2jy  5:  C*(  ||u/  —  vm  II  T  I vl  |2jT  “f"  \vm  I2 ,T  )  > 

we  conclude  that  {u*,}  is  a  Cauchy  sequence  in  J(T)  with  respect  to  the  norm  ||-||2T.  Thus 
there  exists  a  v*  £  J{T)  such  that  linn,-^  Vk  =  v*.  According  to  (2.20),  v*  must  be  a  piecewise 
polynomial  of  degree  not  more  than  1.  By  the  embedding  theorem  of  the  Sobolev  space,  v*  also 
satisfies  the  last  two  equations  in  (2.16),  and  from  the  proof  of  Lemma  2.2,  we  have 


(/rvu*--/?+vu*+).ni^ 


(rvr-rvP+)-nA 

J  DE 

-  1  f  (P~Vv*~  -  /3+Vu*+)  •  nrds  =  0, 

JDEds  JrnT 


Hence  we  have  v*  £  Sh(T)  with 


v*(A)  =v*(B)  =v*(C)  =0. 


By  Lemma  2.1,  v*  =  0  which  contradicts  with  (2.20),  and  this  implies  that  (2.19)  should  be  valid. 
Both  (2.18)  and  (2.19)  plus  some  simple  calculations  show  that  |||  •  |||2,t  is  a  norm  equivalent 

to  ||  •  ||2,T- 


For  any  u  £  J(T),  we  let  Inu  £  Sh{T)  be  its  interpolant  in  Sh(T)  defined  by 
u(A)  =  Ihu(A),  u(B)  =  Ihu(B),  u(C)  =  Ihu(C), 

To  investigate  the  accuracy  of  IhU,  we  will  use  the  reference  triangle  T  whose  vertices  are 


If  the  interface  triangle  T  has  the  following  vertices 


then  T  is  the  image  of  the  T  under  the  affine  mapping  below: 

F(x,y)  =  *(})  +  (£),  (*,£)*  £T, 
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where  the  matrix  M  is  defined  as 

M  _  (  X 2-X1  X3  ~  X!  \ 

\  2/2  -  2/i  2/3  —  /  ' 

Figure  3  gives  a  sketch  of  T  and  its  image  T  under  the  affine  mapping  F(x,y).  For  each  function 
v(x,y)  defined  on  T,  we  can  have  a  function  v(x,y)  defined  on  T  by  letting  v  =  v  o  F.  The 
definition  of  Sh{T)  above  can  be  extended  to  define  Sh(T).  The  following  lemma  states  that 
the  standard  norm  relationship  between  v  and  v  hold  even  when  T  is  an  interface  triangle. 


A  C  =  (0,l)t 


B  C  A  =  (0,0)*  E  -B  =  (1,0)* 


Figure  3:  The  reference  element  T  and  its  image  T  under  the  affine  mapping  F{x.  y). 


Lemma  2.4  For  a  given  (either  interface  or  non-interface)  element  T  and  non-negative  integer 
k  <  2,  there  exists  a  constant  C  such  that 

\v\kj,  <  C  || 1| k  |det  M\1/'2  \v\kf, 

\v\ktf  <C\\M\\k  |det  M\~1/2  \v\k/r 

Proof.  The  results  follow  from  arguments  similar  to  those  used  in  proving  the  same  inequalities 
for  functions  in  standard  Sovolev  spaces,  see  for  example  [9,  16]. 


The  following  theorem  describes  the  approximation  capability  of  Sh{T). 

Theorem  2.1  Let  T  be  an  interface  triangle,  then  for  any  u  £  J(T)  we  have 
(2.21)  \\u  —  Ihu\\m  T  <  Ch2~m  \u\2  T  ,  0  <  m  <  2, 

where  h  is  the  length  of  the  longest  edge  ofT. 

Proof.  First,  we  consider  the  case  with  the  reference  triangle  T.  By  Lemma  2.3, 
u-Ihu  .  <  \\\u-Ihu\\Lf 

2 ,T 

=  C{u-Ihu  ;  +  u(A)-ihu(A )  +  u(B)  —  Ihu(B)  +  u(C)  -  Ihu(C)  ) 

2,T 

=  C  u  —  Ihu  =  C  luL  . 

2  ,T 
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Then  we  can  obtain  (2.21)  by  Lemma  2.4  and  the  standard  homogeneity  argument  [6,  16]. 


Now  we  use  the  partition  Th  to  define  an  immersed  finite  element  (IFE)  space  5^(0).  We 
first  define  a  nodal  basis  function  <j>(x,y )  for  each  node  (xnjVnY  of  Th  piecewisely  such  that 
^(xnjVn)  =  1  but  zero  at  other  nodes,  and  (f>\ t  £  Sh(T)  for  any  triangle  T  £  Th-  Here  Sh{T) 
is  the  usual  space  of  linear  polynomials  when  T  is  a  non-interface  triangle,  or  the  immersed 
finite  element  space  on  T  introduced  in  Section  2.2  when  T  is  an  interface  triangle.  Then  we 
define  Sh(fi)  as  the  span  of  these  nodal  basis  functions,  and  it  is  easy  to  see  that  Sh{ H)  has  the 
following  properties: 

•  For  a  partition  Th ,  the  IFE  space  S^O)  has  the  same  number  of  nodal  basis  functions  as 
that  formed  by  the  usual  linear  polynomials. 

•  For  a  partition  Th  fine  enough,  most  of  its  triangles  are  non-interface  triangles,  and  most 
of  the  nodal  basis  functions  of  the  IFE  space  Sh(ty  are  just  the  usual  linear  nodal  basis 
functions  except  for  few  nodes  in  the  vicinity  of  the  interface  T. 

•  For  any  cj)  £  S/^fi),  we  have 

(2.22)  0|n\n-  £ 

where  O'  is  the  union  of  interface  triangles. 

Now  we  consider  a  function  u  satisfying 

(2.23)  u  £  C(H),  u\QS  £  H2{ns),  s  =  ~,+ 
and 

(2.24)  (/rVu“ -/?+Vu+)  •  n  =  0, 

on  T,  and  we  define  its  interpolant  IhU  in  the  IFE  space  S/,(fi)  by 
Ihu(x,  y)  =  u(x,  y),  if  (x,  y)  is  a  node  of  Th- 
From  Theorem  2.1  we  can  easily  obtain  the  following  error  estimate  for  Tu. 

Theorem  2.2  Assume  that  u  satisfies  the  conditions  (2.23)  and  (2.24),  then 

(2.25)  \\u  —  Tu ||q  q  +  h\\u  —  Ihu\\itQth  —  Ch  ||u||2|^  , 
where 


II  m,Q,h 


y:  iMim.r  • 

TeTh 


3  A  finite  element  solution  and  its  error  estimates 

It  is  obvious  that  the  finite  element  space  Sh{ty  introduced  in  the  last  section  is  not  in  the  space 
to  which  the  solution  of  the  interface  problem  belongs.  A  function  <f>  of  5/,(0)  is  continuous  in 
the  set  of  the  union  of  non-interface  triangles.  The  possible  discontinuous  points  of  <j>  are  on  the 
interface  T  or  on  edges  of  an  interface  triangle  which  intersect  with  T  somewhere  between  its 
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vertices.  Therefore  the  finite  element  method  based  on  Sh{ 0)  is  nonconforming.  To  introduce 
the  finite  element  method  we  consider  the  following  bilinear  form: 

(3.26)  ah(u,v)  :=  /  /3Tu  •  Vvdxdy,  for  all  u,  v  G  H1^)  © 

TeTh  ^TeTh 

and  let 

(3.27)  ||u||ft  :=  ^/ah(u,u),  for  all  u  G  Hl{£l)  ©  5/,( fi), 

which  is  a  common  quantity  used  in  the  error  estimation  of  nonconforming  finite  elements,  see 
for  example  [4].  Also  it  is  obvious  that 

IML  ^  C  IMIi,n,/>  • 

By  simple  calculations  we  can  show  that  this  bilinear  form  has  the  usual  boundedness  and 
coercivity  as  described  in  the  lemma  below. 

Lemma  3.1  There  exists  a  positive  constant  C  such  that  for  any  u,v  G  H1^)  ©  5ft  ( fi),  we 
have 

(3.28)  \ah(u,v)\  <  C\\u\\h\\v\\h, 

(3.29)  ah(u,v)>C\\u\\2h. 

Also,  ||-||ft  is  a  norm  of  the  space  Hq(H)  ©  5/,o(0).  Here  Ho{&)  ©  5fto(0)  consists  of  functions 
in  771(J2)  ©  5/, (fi)  which  are  zero  on  9fi. 


For  the  interface  problem,  we  now  define  its  immersed  finite  element  solution  as  a  function 
uh  G  5ft0  satisfying 

(3.30)  ah(uh,vh)  =  (/>«/>)>  Vvh  G  Sho- 


To  aid  the  error  estimation  for  uh,  we  let  Th [  be  the  collection  of  all  the  interface  triangles,  and 
let  E'h  be  set  of  all  edges  in  Th  which  meet  the  interface  T  between  their  vertices.  Each  element 
T  G  Tl  is  separated  into  two  pieces  T~  and  T+  by  the  interface  T  such  that  T~  c  0“  and 
T+  c  fi+. 

We  first  consider  the  consistency  error  in  the  weak  form: 

(3.31)  Lu(vh)  :=  ah(u,vh)  -  (^f,vh^  . 

Here  u  is  the  solution  of  the  boundary  value  problem,  and  Vh  is  a  function  in  Sho-  From  the 
Green’s  formula,  we  have 


(3.32)  Lu  (vh) 

(3.33) 


E  E  l 

■,e£'h  e  t eT' 

e  (jy£^d‘+iT 


[P^vh]ds, 
Tnr  on 


du 

dn 


vh]ds  + 


\'BtnVh\dS: 


where  ^  is  the  normal  derivative.  For  the  terms  in  the  first  summation,  we  assume  that  e  is 
the  common  edge  of  two  interface  triangles  Tl  and  Tr,  see  the  sketch  in  Figure  4.  Then  we  have 

<334)  /TsH*  =  l- 
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Figure  4:  Two  adjacent  interface  triangle  whose  common  edge  is  cut  into  two  pieces  e  and  e+  by 
the  interface  T. 


because  has  no  jump  in  ( Tl  U  Tr)  n  f 1  .  Similarly, 
(3.35)  =  /+  ~ 

Here  vlh  —  Vh  |t;  ,  i  =  h  r ■  For  the  second  term,  we  have 


(3.36) 


-v+)ds. 


With  these  preparations  we  can  now  derive  an  estimate  for  the  consistency  error  in  the  following 
theorem. 


Theorem  3.1  Assume  that  the  solution  u  of  the  interface  problem  is  in  C(Q)  n  H2(  12).  Then 


(3.37) 


ah(u,  vh)  -  ( f ,  vh ^  <  Ch  \u\2  U  \\vh\\h  , 


for  any  vh  e  Sh(Q). 

Proof.  According  to  the  discussion  above,  we  need  to  estimate  the  quantities  in  (3.34)  -  (3.36). 
For  (3.34),  we  note  that 


Jh(A)=vrh(A) 

where  A  is  the  common  vertex  of  T1  and  Tr  in  fl-.  Then 


(3.38) 


= 

J 

f  ~Vlh(A)  +vrh{A)  -vrh)ds 

< 

J 

[j^«-vlh(A))ds 

+ 

L* 

vh(A)  -  vrh)ds 
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Let  Xs  =  (x,  y)  be  a  point  on  e  .  The  difference  vlh(Xs)  —  vlh(A)  can  be  written  in  the  following 
integral  term 


dvi 


vlh(Xs)-vlh(A)  =  f_^dr), 


dv[ 


where  is  the  tangential  derivative  of  vlh  along  the  line  AXS.  Therefore,  by  the  trace  theorem 
and  the  fact  that  vlh  is  polynomial  of  degree  1, 


J  ~  vlh(A))ds 


f  „du  f  dvlh 

< 

L’tohc-i*** 

<  8 

<  P~hl/2 


[ 

du 

(f 

dvlh 

L 

dn 

\Jax7 

dt 

l.Ml. 


dvi 


di]  J  ds 
2  \  1/2 


<  p-h1'2  (l 


dvi 


dt 
2  \  1/2 
ds  | 


h  '  ds  J  ds 


J 


<  B  h 


dt 


L 


du 


dn 


ds 


<  cp~h\\Wh\\ 

i,T‘nn-  l|V«|| 

l-.T'nn- 

—  Cfi  h  Iul2,7^nn-  ■ 


Similarly, 

and  we  have 
(3.39) 


j  P^(vrh(A)  ~  vh)dt 


<  Cf3  h\vh\1T,,nQ-  |w|2)x' 


nn-  ’ 


f  dn,Vh^S  —  CP  M  M2,T'nn-  +  luftli,Trnn-  \u\2,Trnn-  )• 


By  the  same  argument,  we  can  have 


(3.40) 


J '+^ dnVh^S  ~  l^li,Tinfi+  \u\2,Tlnn+  +  I i,r^nr2+  M2,T*nQ+  )■ 


We  now  turn  to  the  term  on  the  interface.  Assume  that  the  interface  T  meets  the  edge  of  T  at 
a  point  D.  Then  by  the  definition  of  the  functions  in  Sh{T)  we  have 


~(D)=v+(D). 


Therefore  we  have 


ITJ,^n]d‘  =  -^D)+<m-,i)d, 

,Tnrrisr<'’» 


=  / 
Jt  r 


+ 


/ 

Jtc\t 


P+^(v+(D)~vt)ds 


=  L  r  isr  {L  Tit11”)  * + JTnr  V*  (/-  a-wd” ) 
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where  XSD  is  the  part  of  T  n  T  between  Xs  =  (x,  y)  and  D.  Repeating  the  same  arguments  as 
those  used  for  (3.39),  we  have 


(3-41)  \L/^]d 

<  Cmax{^“,  P+}h  ^ 


h  l.'ITil  l  l2,Tnn 


-  +  \v 


h  li,Tnn+  l“l2,Tnn+ 


)■ 


Finally,  the  estimate  of  (3.37)  is  obtained  by  putting  estimates  (3.39),  (3.40)  and  (3.42)  in  (3.33). 


Now,  we  are  ready  to  derive  an  error  estimate  in  the  H-H^  norm. 

Theorem  3.2  Assume  that  the  solution  u  of  the  interface  problem  is  in  (7(0)  fl  FP(O).  Then 
the  error  of  the  immersed  finite  element  solution  uh  has  the  following  estimate  for  a  constant 
C  >  0: 


(3.42)  \\u-uh\\h  <  Ch\u\2>a. 

Proof.  From  the  basic  Lemma  of  Strang  [4],  we  have 


||w  -  uh ||  <  <7 


(  .  t  I,  n  ,  l^«(^ft)|N\ 

mf  sup  -r - n —  . 

^VhGSho(^)  Wh€Sho(ty  II^hI/i  J 


Then  the  estimate  of  (3.42)  is  obtained  by  applying  Theorem  2.2  and  3.1  to  the  above. 


We  now  turn  to  the  error  estimation  in  L2  norm.  Note  that  for  any  g  £  L2( 0),  the  auxiliary 
problem 

(3.43)  /  /3V0  ■  Xvdxdy  =  /  g  v  dxdy,  for  any  v  e  H1( 0) 

Ja  Jn 

has  a  unique  solution  in  (7(0)  fl  H2( 0)  satisfying 
(3-44)  Ua\kn<C\\g\\0,Q- 

We  let  (j>h  be  the  immersed  finite  element  solution  of  this  auxiliary  boundary  value  problem. 
Then,  from  Theorem  3.2  we  have 

(3.45)  1109  -  0/>  1 1  <  Ch  |09  |2)Q  <  Ch  y|on  . 

In  addition,  the  estimate  in  Theorem  3.1  is  still  true  if  we  substitute  Vh  by  a  function  in 
((7(0)  n  H2(n))  ©  Sho-  Thus 


(Lhiy*  'U'hi  ) 

(u  -  uh’9^o 

= 

Oft (09,  U  -  uh)  -  (g,u-  Uh)o 

< 

Ch  |09  l2,fi  llM  —  Uh\\h 

(3.46) 

< 

Ch2  ||5||0)n  \u\2^  . 

Q'hi.'U'i  4*g  0/i)  ^ 

f )  09  —  0ft^ 

< 

Ch\u\2Q\\<f>g  -  KWh 

(3.47) 

< 

Ch  |w|2)n  1109 ll2,o  —  Ch  |w|2)q  ||3||0)n 

These  preparations  lead  to  an  error  estimate  in  the  L 2  norm  in  the  following  theorem. 
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Theorem  3.3  Assume  that  the  conditions  of  Theorem  3.2  are  satisfied.  Then  the  error  of  the 
immersed  finite  element  solution  uu  has  the  following  estimate  for  a  constant  C  >  0: 

(3.48)  \\u-uh\\0U<Ch2\u\2U. 

Proof.  By  the  generalize  Aubin-Nitsche  lemma  [4],  we  have 


llM  -«fc||0,n  <  —  {  \\u~uh\\h  || 4>g  -  4>h\\h  +  \ah{u  -uh,4>g)  -  (u-uh,g^ 

+  Gift  (it,  (j)g  ~  4>h)  ~  ^/,  (f>g  ~ 


The  estimate  (3.48)  is  then  obtained  by  applying  Theorem  3.2,  (3.45),  (3.46),  and  (3.47)  to  the 
above. 


4  A  modified  scheme 

The  immersed  finite  element  space  introduced  in  Section  2  involves  piecewise  polynomials  defined 
over  subsets  with  curve  boundaries  in  interface  triangles.  While  this  does  not  cause  any  difficulty 
in  some  of  its  applications  such  as  the  finite  volume  element  method  [10],  we  would  like  to 
have  a  finite  element  formulation  in  which  the  integrations  are  carried  out  over  approximate 
triangle/quadrilaterals. 


B  C  B  C 


Figure  5:  Subsets  T  and  T+  formed  by  the  interface  T  in  an  interface  triangle  are  modified  to 
T*  and  T*+ 

Let  us  consider  a  typical  interface  triangle  T  6  Th  such  that 

t  =  t~  u  t+,  t~=t  nr,  T+  =  Tnfi+. 
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As  before,  we  use  A,B,C  to  denote  the  three  vertices  of  T  with  A  6  0  ,  B,C  £  fl+,  and 
assume  that  the  interface  T  meets  two  edges  of  T  passing  the  vertex  A  at  D  and  E.  Then,  T~ 

is  a  triangle  with  a  curved  edge  ED ,  and  T+  is  a  quadrilateral  with  the  same  curved  edge  as 
T~ .  On  the  other  hand,  this  triangle  can  also  be  separated  into  T*~  and  T*+  by  the  straight 
line  between  D  and  E,  see  the  sketch  in  Figure  5.  In  the  immersed  finite  element  space  S/,(0), 
we  introduce  a  new  bilinear  form 


a*h(vh,wh)  :=  \ 

r. rtz'T,  \  ti  •'T 


/3Vu/j  •  Vwudxdy 


(4.49) 


TeTh\B 

Ter: 


T— 


/3  VuJ)  ■  VtoJ  dxdy  +  /  /?+Vu£+  •  Vw*h+dxdy 

Jt*+ 


and  a  new  linear  functional 


fZ(wh)  :=  'Yj  fwhdxdy+  ^  [  f 
TeTh\B  reTl 


T*  ~ 


fw*h  dxdy  +  /  fw*h+ dxdy 

JT*+ 


where  for  any  function  Wh  €  Sh(fi),  we  let  w £  be  a  piecewise  linear  function  such  that 


and 

(4.50) 


wh\r(x,y) 


-{ 


*+ 


~(x,y), 

{x,y), 


C x,y )  e  T*  , 
(*,!/)  GT*+, 


wl  (z,2/)  =  «v,  (x,y)tT*  nr  , 
w*+(x,y)=w+,  (x,y)&T*+nT+. 


Then  we  define  the  modified  immersed  finite  element  solution  of  the  interface  problem  as  a 
function  u*h  £  Sho  (H)  satisfying 


(4.51)  a*h(u*h,wh)  =  f*(wh),Vwh  £  Sh0(Tl). 


It  is  obvious  that  a*h(vh,  Wu)  is  a  continuous  bilinear  form  in  Sh{D)  with  respect  to  the  norm 
||-||ft.  For  further  discusssions,  we  use  T  to  denote  the  subset  in  an  interface  triangle  T  £  Th 
formed  by  the  curve  T  fl  T  and  the  line  DE,  and  let 


(4.52) 


T*~  =  T  n  T*~, 

rjn*-\-  _  ^  f)  + 


see  the  sketch  in  Figure  6.  Then,  we  can  show  that  this  bilinear  form  has  the  usual  coerciveness 
as  stated  in  the  lemma  below. 


Lemma  4.1  Assume 


max  ■ 


f*+ 


rp*- 


|T-|  ’  |T+| 


<C,  *  — 


for  a  constant  0  <  C  <  1 .  Then,  there  exists  a  positive  constant  C  such 


a*h(vh,vh)  >  C\\vh\\l 


for  any  vh  £  Sh(Tl). 
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A 


Figure  6:  Sub-sets  formed  by  the  line  DE  and  the  interface  between  D  and  E. 


Proof.  For  each  interface  triangle  T,  we  have 


/ 

JT*~ 


•  Vv!  dxdy  + 


-  L 


8  Vv*h  ■  Vv*h  dxdy  + 


Vv*h+dxdy 


/ 

Jt * 


/?  dxdy 


> 


> 


[  p+Vv 

Jt*+ 

/  P~Vv*h~  ■  Vv*h~dxdy  - 

1  — 

+  [  P+Vv*+ -Vvl+dxdy-  [  P+Vv*+  ■  Vv*+dxdy  +  [  P+Vv*+  •  V« 

Jt+  Jf •-  Jf*  + 

/  P~Vv^  ■  Vv^dxdy  +  /  ,d+V^-V^dxdy 
Jt-  Jt+ 

-  /  P~ Vvf  ■  Vvfdxdy  -  /  /3+Vv^"  •  Va^"dxdy 

Jt*+  Jt*- 

/  P~ ■  Vv^dxdy  +  /  /3+V^-V^dxdy 
Jt~  Jt+ 

-  C  f  P~Vv~  ■  Vv~dxdy  -C  f  /3+Vv+  ■  Vv+dxdy 

Jt-  Jt+ 


dxdy 


=  (1  -C) 
Therefore, 


/  8  Vvh  -  Vvhdxdy  +  /  P+Vv^  ■  Vv^dxdy 

Jr  -  lr+ 

a*h(vh,vh )  >  (1  -  C)ah(vh,vh)  =  C 


2 

/l  ’ 


The  following  lemmas  are  needed  to  derive  an  error  estimate  of  the  modified  immersed  finite 
element  solution  u*h  €  Sho- 

Lemma  4.2  There  exists  a  constant  C  such  that  for  any  Wh  €  Sho{ty  we  have 
(4-53)  I fh(wh)  -  AK)|  <  Ch\\f\\0>n  IKH,,. 
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Proof.  For  each  interface  element  T  e  7^, 

/  fw*h~dxdy+  /  fw*h+dxdy  =  /  fwfdxdy  +  /  fw*h~  dxdy  -  /  fwfdxdy 

Jt*-  Jt*+  Jt-  Jt *-  7t*+ 

+  /  fw^dxdy-  /  fw^dxdy+  /  fw*h+dxdy 
Jt+  Jt—  Jt*+ 

=  /  fwfdxdy  +  /  +  /  f  (w*h~  -  w£)  dxdy 

Jt-  Jt+  Jt*~ 

[_  f{w*h+ ~wh)dxdy. 

Jt*+ 


+ 


Hence,  by  the  definitions  of  fh(wh)  and  f^{wh)  we  have 


(4-54)/h(«;h)  -  ^  f  /K  ~wt)dxdy+  f  f  (w*h+  -  wh)dxdy 

TeV  Ut—  Jt*+ 


Take  a  point  (x„,  yn)  on  the  line  D.E  such  that  the  line  connecting  (x,  y)  €  T*  and  ( xn ,  yn )  is 
perpendicular  to  DE.  Then 

dw~^~  dw~^~ 

wi(x,y)  =  w^(xn,yn)  + -^(xn,yn)(x  -  xn)  + -^-(xn,yn)(y  -  yn) 

=  W ^  ( Xn ,  y„)  +  ~Q^(Xn,  Vn)(x  -  Xn )  +  -^-(x„,  Vn)(y  ~  Vn), 


and 


Similarly, 


[  f(w*h  -w+)dxdy  <  hf  |/|  ( 
Jt*-  Jt*~  \ 

f  i/i 

Jt+ 


<  h 

<  Ch 


'9w+  dwj 
dx  dy 

dwt  ,  dwt  \ 

dx  dy  J 


0 ,T+  \\Wh  \\h,T+  ■ 


dxdy 

dxdy 


f.  fiw*h+  ~wh)dxdy 
Jt*+ 


—  Ch  ||/|lo,T-  Wh  h,T~ 


The  result  of  this  lemma  is  then  obtained  by  applying  these  estimates  in  (4.54). 


Lemma  4.3  There  exists  a  constant  C  such  that  for  any  Vh,Wh  of  Sho{ty,  we  have 
\a*h(vh,wh)  -  ah(vh,wh) I  <  Ch  IIv/jII^  llwfcllfc  ■ 

Proof.  Again,  first  consider  a  typical  interface  triangle  T  £Tf,  and  we  have 


[  /T V<“  •  Vw*~  dxdy  +  [  /3+Vv*+  ■  Vw*+dxdy 

Jt *-  Jt*+ 


-L 


/ 

Jt*— 


ft  Vvh  -  Vwh  dxdy  +  /  ft  \7vh  •  S7wh  dxdy  - 


r  , 

Jt*+ 


P~Vv 


*  VffiJ  dxdy 
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+  f  /?+Vu*+  ■  Vw*+dxdy  -  f  B+Vv;+  ■  Vw*+dxdy  +  [  /?+' Vv*+  ■  Vw\ 

Jt+  Jt *-  Jf*+ 

=  [  P~Vv*h- ■Vw*h~dxdy+  f  /3+Vv*h+ -Vwl+dxdy 

Jt+ 

l 


dxdy 


L 


P  wh  dn  1  h  dn 


ds  + 


fdT*  + 


p  Wh  dn  1  Wh  dn 


ds, 


and  this  gives 

a*h(vh,wh)  =  ah(vh,wh)+  ^  / 

TeT'  ^ 9 


ST*+ 


R-w*~  9v*h  -  B+w*+  dv*h+ 
P  Wh  dn  P  Wh  dn 


ds 


(4.55) 

By  definition, 


+ 


/. 

JdT*  + 


dn 


dn 


ds 


Vh  dn 

on  DE.  Thus,  from  the  proof  of  Lemma  2.2, 


r»r^-«+^= o, 


dn 


j. 

JdT*  + 


dv*~ 

B~w*~ _ - _ B+w*+ 

P  Wh  dn  P  h 


dv 


*+ 


dn 


ds 


=  h 


aT*+nr 


=  [ .  [p  K" 

JdT*+ nr  L 


ds 


.  dv 


*+ 


+B+(w*+(D)-w*+y-gn 


ds. 


where  D  is  a  point  on  dT*+  n  DE.  Following  the  same  argument  as  that  used  for  Theorem  3.1, 
we  have 


dvV 


ds 


B~{w*h~  -  w*h~(D))- 
aT*+nr  un 

~  dv*+ 

B+(w*+-w*+(D))^-ds 
8T*+ nr  un 


,dv~h 


ds 


B  K  -  wh  (D))  dr, 

dT*+n  r  un 

B+(w+-wt(D))-^ds 
aT*+ nr  un 


<  Ch  I  w 


h  ll  ,T~ 


<  Ch  \wt 


h  ll,T+  I  h 


Then  the  result  of  this  lemma  is  obtained  by  applying  these  estimates  in  (4.55). 


Finally,  we  can  derive  an  error  estimate  for  the  modified  immersed  finite  element  solution  u*h  € 
Sh o  in  the  following  theorem. 

Theorem  4.1  Assume  that  the  solution  u  of  the  interface  problem  is  in  C(Q)  PI  and 

that  the  condition  in  Lemma  f.l  is  satisfied.  Then  the  modified  immersed  finite  element  solution 
u*h  has  the  following  estimate  for  a  constant  C  >  0: 

\\U*h  —  u\\h  <CMM2,n  +  ll/llo,fi)- 

Proof.  Let  vh  be  arbitrary  function  in  Sh o,  then,  from  Lemma  4.1,  we  have 

C\\u*h-vhfh  <  a*h(u*h-vh,u*h-vh) 

=  ah{uh  -  vh,u*h  -  vh)  +  {ah(vh,u*h  -  vh)  -  a*h(vh,u*h  -  vh)} 

+  {fh«~vh)  ~  fh(K~vh)}- 


ll  ,T~ 

1  ,T+  ’ 
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Thus. 


C\\u*h-vh\\h<C\\uh-vh\\h+  sup 


\ah(vh,Wh)  ~  a*h(vh,Wh)\ 
Wh^-Sho  IKIL 


+  sup 

wh€Sho 


I  fh(Wh)  ~  fhj'Wh)  | 

Ik*  1 1 /, 


Taking  Vh  =  «ft  in  the  above  and  applying  Lemma  4.2,  4.3  and  Theorem  3.2  lead  to 
C\\u*h-uh\\h  <  Ch{\\uh\\h  +  ||/||0)n) 

<  Ch(h  |rtft|2)Q  +  IM|ft  +  ||/Ho,q)- 

The  estimate  above  and  the  triangle  inequality  lead  to  the  result  of  this  theorem. 


5  Numerical  Examples 


We  now  present  some  numerical  results  to  illustrate  features  of  the  immerse  finite  element 
method.  Because  of  its  simplicity,  we  only  present  results  obtained  by  using  the  modified  IFE 
method  based  Cartesian  partitions  to  solve  the  interface  value  problem  defined  by  (2.5)-(2.8) 
in  the  rectangular  domain  f]  =  (—1, 1)  x  (—1, 1).  The  interface  curve  T  is  a  circle  with  radius 
r0  =  7t/6.28  which  separates  f2  into  two  subdomains  fi-  and  fi+  with 

fi“  =  {(x,2/)  :  x2+y2<rl). 

The  boundary  condition  function  g(x,y)  and  the  source  term  f(x,y)  are  chosen  such  that  for 
a  =  3, 

if  r  <  r0, 
r%,  otherwise, 


(5.56) 


u(x,y)  =  < 


l  J- 


+ 


is  the  exact  solution  of  this  interface  problem  for  a  coefficient  function  j3 (x,y)  such  that 


13  ,  (x,y)  e  n  , 
P+,(x,y)  €  n+. 


Here 


=  \/ x2  +  y2 


the  domain  n  and  the  curve  T  are  sketched  in  Figure  7  together  with  a  typical  partition  for  our 
numerical  results.  We  will  use  the  following  quantities  to  describe  the  errors  in  an  IFE  solution: 

{lbft-«Ho>  when  s  =  0, 

i|uft-u||ft,  when  s  =  1, 

maxj£ivh  |w(xj)  —  w(xj)| ,  when  s  =  oo. 


where  Nh  is  the  set  of  nodal  points  of  the  partition.  These  quantities  represent  the  error  in  the 
usual  L2,  H1 .  and  the  discrete  L°°  (at  the  nodal  points)  norms,  respectively. 

Table  1  contains  actual  errors  of  the  IFE  solutions  with  various  partition  size  h  for  the 
boundary  value  problem  with  the  coefficient  function: 


(x,y)  €  ^  , 
(x,y)  €  ft+. 
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Figure  7:  A  typical  grid  for  the  boundary  value  problem. 


By  simple  calculations,  we  can  see  that  the  data  in  this  table  obey 

II Uh  ~  «Ho  ~  Ch2, 

II uh  ~  u\\h  fa  Ch 


as  predicted  by  the  error  estimates  in  the  previous  sections.  Similar  behavior  of  the  IFE  solutions 
can  be  seen  in  Table  2  for  the  boundary  value  problem  with  the  coefficient  function: 


(x,y)  £  n  , 
(x,y)  €  n+. 


The  IFE  method  also  works  well  for  the  case  in  which  the  coefficient  function  has  a  large  jump, 
see  Table  3.  The  errors  in  this  group  of  computations  are  also  within  the  prediction  of  the  error 
estimates. 

The  data  in  Tables  1  and  2  seem  to  suggest  the  second  order  convergence  of  the  IFE  method 
in  the  discrete  L°°  norm.  However,  results  in  Table  3  indicate  that  this  can  not  be  true  for 
all  the  cases.  The  question  under  what  conditions  the  IFE  solution  can  have  a  second  order 
convergence  in  the  L°°  norm  is  still  open. 


h 

ei  (h) 

e0(h) 

^OO  (^) 

1/8 

1.2761  x  lO^1 

1.0063  x  10~2 

4.6058  x  10~3 

1/16 

6.3613  x  10~2 

2.5239  x  10^3 

1.0407  x  10^3 

1/32 

3.1786  x  10~2 

6.3079  x  10~4 

2.7979  x  10^4 

1/64 

1.5868  x  10~2 

1.5789  x  10~4 

6.4752  x  10~5 

1/128 

7.9317  x  10~3 

3.9451  x  10~5 

1.6799  x  10^5 

Table  1:  Numerical  results  for  the  case  when  (i  =  1,  ,d+  =  2. 
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h 

ei(h) 

e0{h) 

^OO  (^) 

1/8 

1.7544  x  10^ 1 

1.9584  x  10~2 

2.3746  x  IQ-'6 

1/16 

8.7530  x  1(T2 

4.9074  x  10-3 

6.5468  x  10~4 

1/32 

4.3729  x  10~2 

1.2292  x  10-3 

1.8158  x  10~4 

1/64 

2.1851  x  10~2 

3.0733  x  10~4 

4.3849  x  10~5 

1/128 

1.0922  x  10~3 

7.6875  x  10~5 

1.2411  x  10~5 

Table  2:  Numerical  results  for  the  case  when  fi  =  2.  =  1. 


h 

ei(h) 

e0{h) 

^OO  (^) 

1/8 

3.1737  x  10~2 

2.7503  x  10~3 

4.1217  x  10~3 

1/16 

1.6464  x  10~2 

7.3846  x  10~4 

9.5219  x  10^4 

1/32 

8.6090  x  10~3 

1.4499  x  10~3 

5.7889  x  10~4 

1/64 

4.2113  x  10~3 

4.3236  x  10~5 

1.0918  x  10^4 

1/128 

2.1154  x  10~3 

9.3966  x  10-e 

4.2197  x  10~5 

Table  3:  Numerical  results  for  the  case  when  /3  =  l,/3+  =  1000. 
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